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Abstract 

' Several topics in the theory of magnetized turbulence are reviewed with application to star 

' formation and the interstellar medium. The density, pressure, and temperature distribution in 

, a turbulent interstellar medium is described in comparison to a medium dominated by thermal 

04 ' instability. The derivation of the empirical theory for the energy spectrum of hydrodynamic 

^ I and MIfD turbulence is outlined, and comparisons are given to numerical models for the 

U . magnetized case. Next discussed is how density fluctuations modify observations of velocity 

fluctuations in line centroids, coincidentally cancelling the effects of projection smoothing. 
The analytic description of velocity structure functions as a function of the dimensionality 
of the dissipative structures is then covered. Finally the question of whether turbulence can 
support against gravitational collapse is addressed. Turbulence both prevents and promotes 
collapse, but its net effect is to inhibit collapse, so it does not trigger star formation in any 
' significant sense. 

o ■ 

^ 1 Observations 

For more than half a century, evidence for the turbulent nature of the interstellar medium 
O (ISM) has been known (von Weiszacker 119431 [TMTl Miinch .iSSSjl . Armstrong, Rickctt, & 
' Spangler H1995|l compiled observations of the ionized ISM that appear to show that the mag- 
O ^1 netized gas has a fluctuation power spectrum consistent with Kolmogorov H1941fl turbulence 
I ■ over scales ranging from 100 pc to small fractions of an AU. 

^ Other observations pointing to the presence of supersonic, magnetized turbulence in the 

ISM include the presence of superthermal linewidths in the diffuse ionized medium, the warm 
Ci and cold neutral gas, and molecular clouds. Line centroid variations also show correlations 
^ with separation suggestive of turbulence. Magnetic field strengths are measured using the 
• ^ . Zeeman effect or, more recently, the Chandrasekhar-Fermi p953fl method (Heitsch et al. 121)011 
^ ; Ostriker, Stone, & Gammie IBDDll Matthews, Fiege, & Moriarty-Schieven 120021 Crutcher 
et al. 12004(1 . Small scale fluctuations are measured using refractive and diffractive pulsar 
scintillation, and fluctuations in rotation measure and dispersion measure. 

This review briefly discusses the density and temperature structure of turbulent gas, as 
a counterpoint to multi-phase models of the ISM. Analytic descriptions of magnetized tur- 
bulence are then treated, as well as numerical tests of those descriptions. The observational 
consequences of density fluctuations on the line centroid velocity fluctuation spectrum are 
briefly noted. Finally the question of whether magnetized turbulence can support against 
gravitational collapse is reviewed. Topics that are not covered that fall under the general 
topic include magnetized turbulence generated by the magnetorotational instability, dynamos 
and fleld structure in turbulent flow, and diffusion of cosmic rays or tracers such as heavy 
elements. 
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2 Density and Temperature Structure 

Classical theories of the interstellar medium have emphasized the thermal phases produced by 
the structure of the cooling curve as a function of temperature (Field, Goldsmith, & Habing 
[TI)6 9,: for a modern example see Wolfire et al. 1995J. At constant pressure, isolated points of 
stable equilibrium can be found on the pressure-density plane, giving phase structure to the 
ISM. However, turbulent flows generate continuous distributions of pressure (Mac Low et al. 
Er)0 4 Avillez & Breitschwerdt 2004), a point also emphasized by Breitschwerdt and Avillez 
in these proceedings. Although thermal instability still acts in such flows, broad regions on 
the pressure-density plane end up occupied, rather than isolated phases. 



3 Scaling Relations in Magnetized Turbulence 

Our understanding of magnetized turbulence rests on the foundations of the Kolmogorov 
1)1941(1 theory for the energy spectrum of incompressible hydrodynamical turbulence. This 
theory requires two assumptions: that energy input at large scale and energy dissipation 
at small scale balance, and that energy transfer between spatial scales occurs only between 
neighboring scales in the intervening inertial range. Or, to quote Richardson, 

Big whirls have little whirls, which feed on their velocity, 
and little whirls have lesser whirls, and so on to viscosity. 

The energy of incompressible turbulence E oc v'^/2, and the crossing time of a perturbation 
over any scale £ is At = £/vi, where vi is the average velocity difference across structures at 
that scale. The assumptions of constant flow between neighboring scales requires that the 
energy flow through any scale 

i^.oc-oc^ (1) 
must be constant. Thus, the velocity must scale as 

ve oc ^1/^ (2) 

In wavenumber space, Vk oc 

We can derive the energy spectrum in wavenumber space 
E[k) in either ID or 3D. In ID, the energy in some portion of the inertial range 

E = J E{k)dk oc vl. (3) 

Substituting the scaling relation derived from equation^ kE{k) oc fc^^/"^, so E{k) oc k^^/^. 
In 3D, £; = / E{k)(Pk, so k^E{k) oc ^nd E{k) oc 

A theory of incompressible MHD turbulence was proposed by Iroshnikov H1963|l and 
Kraichnan H1965fl based on the idea that energy transfer between scales would still be lo- 
cal, but transferred by interactions between oppositely directed Alfven waves rather than 
vortices. Isotropy was still assumed in their derivation, however, as was a small ratio of fluc- 
tuating to mean magnetic fleld strength. The interaction time between Alfven waves at any 
scale 6t = 1/va is determined by their propagation speed v\ = / (47rp). The change in any 
energy from each interaction 

AE oc -^^t oc veieSt. (4) 
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The perturbation crossing time At determines the change in velocity ve = vi/ At oc vj/£, so 

AE (X ^ — cx — . (5) 

e VA VA 

The time required for energy to pass through an incoherent cascade can be shown to be 



tcoc^—j 5tc.-^, (6) 

so the energy flow E cx v^/tc oc vf/{tvA)- If the energy flow through each scale is constant, 
then Vi oc in this case. The energy spectrum in ID can then be derived analogously to 
the hydrodynamic spectrum, yielding E{k) oc 

However, MHD turbulence in the presence of a mean field with energy of the same order 
as the flow is not isotropic. Goldreich & Sridhar l| 19951 11997|l argued that the anisotropy 
can be taken into account by treating the energy cascade along parallel lines as Alfven wave 
interactions with time scale t\\ — tjvAi following Iroshnikov and Kraichnan, but treating 
the perpendicular cascade with hydrodynamic interactions that merely stir field lines with 
time scale t^ = Ijvi. If the cascades in the two directions maintain what is referred to as 
critical balance and proceed at equal rates, so that t\\ = t±, then k^^VA = kj_Vk- From the 

— 1/3 

hydrodynamic cascade, we know that Ufc oc fcj^ , so the anisotropy of the turbulence can 

2/3 

be derived to be scale dependent: fc|| oc kj^ . The parallel cascade is predicted to still be 
an Alfven cascade with E\\{k) oc fcy^^, while the perpendicular cascade looks more like a 

Kolmogorov cascade with E±{k) oc k^^. The ID spectrum is predicted to be E{k) oc fc^/^. 

Maron & Goldreich H2()()l(l confirmed the anisotropy predicted, but found a ID spectral 
index of 3/2. Cho & Vishniac p?(¥)jl . and Maron, Cowley, & McWiUiams on the other 

hand, do find the predicted 5/3 spectral index. The difference appears to be the strength 
of the initial field. If the fluctuations in the magnetic field are small compared to the mean 
field, the local Alfven cascade dominates, as found by Maron & Goldreich 1)200 l|l . If the 
fluctuations are larger, though, the perpendicular hydrodynamic cascade takes effect, as seen 
by Cho & Vishniac H2000|l and Maron ct al. (j2004|l . Finally, if the fields are weak compared 
to the flow, as in the case of the turbulent dynamo, the magnetic cascade is no longer local, 
and the power spectrum of the field is dominated by high-fc modes, while the velocity power 
spectrum remains dominated by the driving scale (e.g. Maron et al. | 



4 Observations of Velocity Spectra 

One diagnostic of turbulence is the fluctuation in velocity centroid from point to point across 
a map. This is a projection of the 3D velocity fluctuation spectrum into 2D. The 3D spectral 
exponent is O = 1/3 in the incompressible hydrodynamical case (eq. EJ, while for shock- 
dominated Burgers turbulence = 1/2 (Burefers ll974|l . The projection to 2D should reduce 
the value of 8 by 0.5 in the incompressible case (von Hoerner [1951 Brunt et al. I2UU3|I . 
However, the observed value of the 2D velocity centroid fluctuation spectrum is close to 0.5 
(e.g. Brunt & Mac Low 1200^ . as if no projection smoothing had occurred. 

Supersonic turbulence drives density fluctuations. Velocity centroids depend on both 
velocity fluctuations and density fluctuations, as the emissivity at any particular velocity 
along a line of sight depends on the density at that point. These density fluctuations add 
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extra power at small scales, increasing the slope of the velocity centroid fluctuation spectrum 
(Lazarian & Esauivel 1200^1 Brunt & Mac Low 1200^ . By a lucky coincidence, the density 
fluctuations in strongly supersonic turbulence cancel the effects of projection smoothing rather 
accurately. This does not hold for transsonic or subsonic turbulence however, and there is a 
transition region in the mildly supersonic regime where the cancellation is incomplete (Brunt 
& Mac Low l^ITO|l . 

5 Velocity Structure Functions 

Velocity structure functions offer another way of statistically characterizing turbulence that 
has proved tractable to analytic description. The structure function of order p at scale L is 

S,{L) = {\v{x + L)-v{x)\n. (7) 

The structure function of order two can be directly related to the energy spectrum E{k). In 
the inertial range of isotropic turbulence, structure functions scale as a power of the length 
scale Sp{L) oc L^(p). The behavior of the power-law exponents C(p) as a function of the order 
p is argued to depend on the dimension D of dissipative structures in the turbulence, along 
with the scaling exponents for v{L) oc and t{L) oc L/v oc L'^ (She & Leveaue 119941 
DubruUe . 1994,1 . The form of the scaling is 

C(p) = e(i-A)p + c(i-sf«), (8) 

where the codimension C = 3 — D, and S = 1 — A/C. In Kolmogorov turbulence, equation ^ 
shows that 6 = 1/3, so A = 2/3. 

In incompressible hydrodynamical turbulence, the dissipative structures are linear vortex 
tubes with dimension D = 1 (She & Leveque I1994|l . Miiller & Biskamp H2UUU|I showed 
that incompressible MHD behaves as if its dissipative structures have dimension D = 2; 
they suggested that these structures are current sheets. Compressible MHD has also been 
examined. Boldyrev H2002|l and Boldyrev, Nordlund, & Padoan H2002|l suggest that the 
dissipative structures have D = 2 and are primarily shocks. Padoan et al. (|200^ show 
that in super-Alfvenic but subsonic to transsonic turbulence 1 < D < 2 depending on the 
sonic Mach number A4ims- However, Vestuto et al. (|2(J(J1I) raise a note of caution, as they 
find steeper power laws for energy spectra than would be predicted by the velocity structure 
function of order 2 under this formalism. 

6 Gravitational Support 

An important question for understanding star formation is whether magnetized turbulence 
can provide support against gravitational collapse. In order to provide significant support, it 
must act like a pressure in all directions, and it must last longer than a free-fall time. Neither 
of these requirements is clearly met by supersonic turbulence in molecular clouds. 

Furthermore, hydrostatically supported objects do not form easily in supersonic turbulent 
fiows. This is because an isothermal object in hydrostatic equilibrium requires a balancing 
external pressure to confine it. Too high a pressure sends it into collapse, while too low 
a pressure allows it to again expand. Vazquez-Semadeni, et al. H2(JU4fl have demonstrated 
numerically that these two paths are taken by virtually all density fluctuations in a magnetized 
turbulent flow, leaving very few objects that might be subject to ambipolar diffusion and loss 
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of magnetic support over long time periods. Rather, the star formation rate appears to be 
limited by the necessity of assembling supercritical regions by moving gas along field lines 
in order to increase the mass-to-flux ratio, and by the global limitation of collapse by the 
turbulence, when that is effective. 

If turbulence decays in less than a free-fall time, it will not provide significant support 
against gravitational collapse. For several decades, Arons & Max H1975|l was interpreted to 
mean that magnetic fields could substantially reduce the decay rate of turbulence. In the 
late 1990's, three groups demonstrated numerically that even magnetized turbulence decays 
quickly (Mac Low et al. lT^ Stone, Ostriker, & Gammie lT^ Padoan k Nordlund lTU^ . 
Mac Low et al. (|1998|l showed that supersonic, trans- Alfvenic or super- Alfvenic turbulence 
decays as E{t) cx t~^, with a resolution study ranging from 32"^ to 256'^ zones. 

There are a couple of exceptions to this quick decay known, although their astrophysical 
relevance remains to be demonstrated. Biskamp & Miiller (jl999|l used 512^ simulations of 
incompressible MHD to show that E (x t^^ for flows with magnetic helicity H = A ■ B — Q, 
but that E cx t^^^^ for flows with significant helicity. Unbalanced cascades of Alfven waves 
have been shown to decay more slowly by Maron & Goldreich l|2U(Jlfl and Cho, Lazarian, & 
Vishniac (^2002)) . However, to reduce decay significantly, Alfven wave fluxes must be as much 
as an order of magnitude higher in one direction than the opposite. 

Mac Low 1)19991 Ol)(J3|l estimates that the dissipation rate for isothermal, supersonic tur- 
bulence is 

e^~il/2)pv^ULd, (9) 
where Ld is the driving scale. The dissipation time for turbulent kinetic energy 

Trf = e/e ~ L/wrms, (10) 

which is just the crossing time for the turbulent flow across the driving scale (Elmegreen 

mm . 

stone et al. (TW) and Mac Low (TW^ showed that supersonic turbulence decays in less 
than a free-fall time under molecular cloud conditions, regardless of whether it is magnetized 
or unmagnetized. The hydrodynamical result agrees with the high-resolution, transsonic, 
decaying models of Porter & Woodward fl992') and Porter, Pouquet, & Woodward H1994|l . 
Mac Low ( 1999j showed that the formal dissipation time Td = e/e scaled in units of the free 
fall time is 

1 /32\^^^ K K 

Td/Tff = hr Tl - 3-9 . 11 

47r^ V 3 / A^rms Mrms 

where ^ = 0.21/7ris the Kolmogorov energy-dissipation coefficient derived by Mac Low H1999)l . 
■A^rms — tirms/cs is the rms Mach number of the turbulence, and k is the ratio of the driving 
wavelength to the Jeans wavelength Aj. In molecular clouds, A^rms is typically observed to 
be of order 10 or higher. If the ratio k < 1, as is probably required to maintain gravitational 
support (Leorat, Passot, & Pouauet I199()|l . then even strongly magnetized turbulence will 
decay long before the cloud collapses, and not markedly retard the collapse. 

Supersonic turbulence both prevents and promotes collapse. It can provide support against 
collapse on scales larger than the driving scale. Classical theories treat small-scale turbulence 
as an isotropic pressure with effective sound speed Cg.cS = (cs + (■f^^)/3)^^^ (Chandrasekhar 
119511 von Weiszacker 1195 J)l . By this argument, supersonic turbulence increases the Jeans 
mass to Mj,,s = (7r/G)3/2p-i/2c3^ff. However, supersonic fiows drive strong density fluctua- 
tions on scales below the driving scale. Positive density fluctuations in isothermal turbulence 
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reach densities p' = pAd'^^^. In these density enhancements, the Jeans mass decreases to 
Mj^eff = (7r/G)'^/^p'~^/^Cg gjj. If we express everything in terms of velocity, however, 

M,,eff = oc Hi (c, + {v')/3f' oc v\ (12) 

Even though turbulence can locally promote collapse, on average it inhibits collapse when 
compared to the same conditions absent turbulence. Thus, although turbulence may give 
the appearance of locally triggering collapse and star formation, it is globally preventing star 
formation, acting as a counterbalance to gravitational instability. Turbulence controls star 
formation in opposition to gravity, not in alliance with gravity. 

I thank J. Maron & A. Schekochihin for useful discussions, the conference organizers 
for partial support of my attendance, and J. Pulte for hospitality during the writing of the 
proceedings. Further support was provided by NSF grant AST99-85392 and NASA grant 
NAG5-10103. Extensive use was made of the NASA Astrophysics Data System in preparation 
of this work. 
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